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Abstract 

We calculate the single-spin and double-spin asymmetry differential cross 
sections for the polarized hadron scattering PP — ► l + P + jet up to 0(a s ) 
by the helicity amplitude method. Numerical results of the differential cross 
sections, which can be used to probe the spin contents of the proton, are 
obtained from several sets of polarized parton distribution functions. 
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I. Introduction 

In recent years, the determination of the spin contents of the nucleon has 
been the subject of a large number of work, which was initiated by the EMC 
experimental results [1] that the spin of the proton does not arise from its 
constituent quarks as expected in the naive quark model; the valence quarks 
contribute very little to the spin of the nucleon. 

Presently, there are three possibilities for explaining the EMC and other 
subsequent experimental results [2,3,4], namely (i) large sea quark polariza- 
tion [5]; (ii) large gluon polarization [6]; and (hi) mediately large sea and 
gluon polarizations [7]. Since a different explanation leads to a different 
set of polarized parton distribution functions, the determination of the po- 
larized parton densities by other experiments like polarized hadron-hadron 
scattering is very important for understanding of the spin contents of the nu- 
cleon. From the experiments, we can study the single-spin and double-spin 
asymmetries. For both initial protons longitudinally polarized, we study 
the double-spin asymmetries. Here the inclusive cross sections for the in- 
coming hadrons' longitudinal polarizations are either parallel or antiparallel. 
From the experimental point of view, it is simpler to consider the single-spin 
asymmetries, where one of the initial hadron beams is longitudinally polar- 
ized and has positive or negative helicity. A nonzero single-spin asymmetry 
implies that some of the parton-parton scatterings involve parity-violating 
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weak interactions. Therefore, the single-spin asymmetries can also be used 
to probe parity- violating parton-parton subprocesses [8]. 

In this paper, we shall study the polarized PP — > l + l~ + jet scattering 
by the helicity amplitude method [9]. Both the single-spin and double- 
spin asymmetries are studied. At high energies, it is expected that virtual 
photon and virtual Z boson equally make important contributions to the 
lepton pair production. We organize the paper as follows. In section II, we 
shall use the helicity amplitude method to calculate the helicity amplitudes 
for the differential cross sections for the polarized hadron scattering. In 
section III, we calculate numerically the single and double asymmetries for 
the scattering. Finally, our conclusion is given in section IV. 

II. Calculation of the Helicity Amplitudes 

In this section, we shall calculate the amplitudes for the single-spin and 
double-spin asymmetry differential cross sections for the polarized PP — ► 
l + l~ + jet scattering up to 0{a s ) by the helicity amplitude method. We 
need to consider the following two parton subprocesses and their charge 
conjugates contributing to the dilepton production 



G + q -> l + l~+q, 



(1) 



q + q->l+r + G. 



(2) 
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We denote M(Ai, A2; A3, A4, A5) as the helicity amplitude, where Ai and A2 
are the helicities of the initial partons, while A3 A4 and A5 are the helicities of 
the lepton pair l + , l~ , and the final parton, respectively. The nonvanishing 
tree-level amplitudes read 

S/f /x \ \ \ x\ / M(Ai, A 2 ;A 3 ,-A 3 , A 2 ), foiG + q^ l + l~ + q , 
M(A 1 ,A 2 ,A 3 ,A 4 ,A 5 ) = | M(Al; _ Ai;A3; _ A3)A5)) iorq + g^i + i- + G . 

The helicity amplitudes for figures l.(a) and l.(b) are M(Ai, A 2 ; A3, —A3, A 2 ), 
whereas those for figures l.(c) and l.(d) are M(Ai, — Ai; A 2 , — A 2 , A3). 

To evaluate the amplitudes, we first define some useful formulae and our 
convention. We denote the positive and negative helicity states by | A±), 
which have the following properties: 

(1±75)|Af> = 0, 

\A + ) C = \A.), 

(A T \B ± ) = -(B T \A ± ), 

(A+ I 7^ I B+) = <B_ I 7 „ I A_> , 

and, by Fierz rearrangement theorem, 

(A+ I 7^ I £ + )<C_ I 7 M I D.) = 2{A + I £>_)<C_ \ B + ) , 

(A. I B + )(C_ I £>+) = <A_ I D + )(C_ I B+) + (A. \ C + ){B_ \ D+) . 



(3) 



(4) 
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The massless spinors with momentum p and helicity A = ± are u±(p), 
v ± (p) > u± (p) an d v± (p) > which satisfy 

p\i(p) = $v{p) = u(p)tf = v(p)$, p 2 = , 
(l±7 5 )t;± = (1=F7 5 )«± = u±(l±7 5 ) = u±(lT75) = 0, (5) 
«±(p)7/u«±(?>) = t>±(p)7/i«±(p) = 2p^. 

For massless momentum p, we can write ^ =| p + )(p + | + | p-)(p- | and use 
the convention 

u±(p) = v T {p) =\p±), 

u±(p) = v T {p) = (p± I , (6) 

Ipt) = \p±Y- 

For simplicity, we shall write 

(p- \q+) = (pq) = -(qp) , 

(q+ \p-) = {pq)* = - (qp)* , (7) 

I (pq) I 2 = 2p-q. 
Finally, the gluon helicities are chosen to be 

4 = j^rzr[\P?M\ + \pt)<pt\] ^ G + q^l+r+q, (8) 
\Pb I Pi I 

4 = , tf±A \Ps)(Pi\ + \pt)(pftt ^ q + q^l+r + G. (9) 
\Pi I P5 I 
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At high energies, both the virtual photon and Z boson contribute to the 
lepton pair production. The helicity amplitudes are listed in Appendix A. 
From the amplitudes, we can obtain the polarized parton differential cross 
sections as follows: 



da^±da^ oc (| M(++; + -+) | 2 + | M(++;- + +) | 2 

± | M(+-; + - -) | 2 ± | M(+-; - + -) | 2 ) fig .ia,ib 
±(| M(+-; + - +) | 2 + | M(+-; - + +) | 2 
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To write the differential cross sections explicitly, we consider the center- 
of-mass frame of the incoming hadrons with momenta 



Pi = 1 o, o, i) , p 2 = ^(l, o, o, -l) 



(10) 



Let the momentum of the virtual photon or Z boson be q. Then the Man- 
delstam invariants are: 

S = (P1+P2) 2 , 

T = {q-P2?, (11) 
U = {q-Pif. 

At the parton-level scattering processes, the parton momenta have the ex- 
plicit forms: 

Pi = X1P1 , 
P2 = X2P2 , 

= ri+pZ-Q", (12) 

Q M = P$ + Pl, 
<f = Pl-Pl, 

where p$ is the momentum of the outgoing parton. The Mandelstam invari- 
ants are: 

s = (p! + p 2 ) 2 = xix 2 S , 

i = (p 5 - P2? = (Q ~ Pi? , (13) 

U = (p 5 -Pi? = (Q -P2) 2 ■ 

In the center-of-mass frame of incoming partons, the absolute squares of 
the sum of the helicity amplitudes for the subprocesses G\ l (pi)+ q\ 2 {P2)-^ 



l \ 3 (P3) l \ 4 (P*)+ 9A 5 (P5) are (cf. Appendix B): 



|M(++; + -+) \i g 

.la,lb 



+ 



8<? 2 e 4 Q 2 
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s s 

IT (Q 2 t — Su) , ^9„- ^0-9 
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(14) 
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|M(+-; + - -) |i g . la , 



lb 



- Mf)2 + M|r|] 
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(17) 



While for the subprocesses q Xl (Pi) + <?A 2 (P2) ^(Pa^G^) + G As (p 5 ), the 
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absolute squares of the sum of the helicity amplitudes read 



v,2 ( 2g 2 g%L 2 q L 2 e 8g 2 e A Q 2 q 

| M( + -; + -+) |fig.lc,ld - r/ 2 _ A#2 A2 T A/ f2 r2 1 + 



+ 
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3 (s-Q 2 ) 2 3 (s-Q 2 ) 2J 



I Af (+-; - + -) || g . lc , ld Re = Le | M(+-; + - -) | 2 g . 1Cjld . (21) 



III. Numerical Results 

In this section, the differential cross sections for the polarized processes 
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will be calculated numerically. For the double polarized hadron process, we 
have 



]T J dxi dx 2 Mi{xi)mx2)[s^^] , (22) 

- — ij 

where y is the rapidity, 

A/i(s) = /<+>(*) -/<">(*), (23) 
is the polarized structure function, and 

dAa, tl = ^(da^ +) -da^). (24) 
For the single polarized hadron process, we have 

(25) 

= ^(< +) + <- } -< + ) -<-)). (26) 

In the rest frame of 7*(^°) [10], the momenta of the lepton pair and 
outgoing parton are: 

P3 = ^(E' — q cosa,q sinO — qsmacosf3cos9 — E' cosas'mO, 
—q sin a sin /3, q cos — E' cos a cos + q sin a cos /? sin 9) , 

P4 = -{E 1 + (/ cos a, g' sin 6* + (/sin a cos /? cos + E 1 cos asin#, (27) 
<7 sin a sin /3, q' cos + E' cos a cos 6> — g sin a cos /? sin 9) , 
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where E' = and q' = s % . The above momenta are in the form 

x M = (t,x,y,z). The momenta of the incident partons thus define the z axis 
and the direction of the j*(Z°) defines the x-z plane. The angles a and 
(5 describe the decay of the 7*(Z°) relative to its axes. The phase space 
integration becomes 

1 d?pz<fip A d?p b A 

jr^nr 5 (P1+P2-Q-P5) 



(2vr)5 2pl <lp\ 2p% 



1 1 ,^ r,~ - ~ ,dQ 2 dtdu . 

: d07r(5(s + t + u - Q 2 )(-^ ) .(28) 



(2vr) 5 16 

To calculate the differential cross sections at the hadron level, we write the 
(double polarized) differential cross section as [11] 

<IA(T - £ /\/.r l r/.r. i A/i.r l iA/S.r. i );.s--' /ArT '' 



dQ 2 dydql J dQ 2 didu 



j dx±dx25(s + i + u — Q 2 )f(xi, x-i) 
dxif(x 1 ,x* 2 ) 
"r+& xiS - ^(Q 2 + g£)5e» 
+ f 1 dx 2 /(xt,x 2 ) ^ 

J Vr+ e ~ v X2S — y/s(Q 2 + (7j.)5e~2/ ' 



where 



and 



- W T + f 5 T = T' (30) 



■^1 — ^ i i l"-U 

x 2 s-Vs(Q 2 + q 2 )2 e y 
x * = x 1 Vs(Q 2 + g 2 )hy-Q 2 
2 Xl s-Vs(Q 2 + q 2 )h-y' 
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Finally, we evaluate the single-spin and double-spin asymmetries defined 
respectively as: 

da ( + ) _ da (-) 
Al ~ da ( + ) + da (-) ' ^ 
_ d£7 (+-) 

= *,(++> + *,(+-> • (M) 
To evaluate A^^, we need polarized parton distribution functions, which 
we shall use those given by Cheng, et al. (CLW) [12] and those by Gliick, et 
al. (GRSV) [13]. To obtain Al, we need to use polarized parton distribution 
functions together with the corresponding unpolarized ones. Specifically, the 
CLW sets should be used with the MRS(A') sets [14], while the GRSV sets 
should be used with the unpolarized GRV sets [15]. 

In the CLW sets, two different cases of polarized parton distributions 
in the gauge-invariant factorization scheme are used: (I) sea polarization 
As(x) / 0, gluon polarization AG(x) = 0; and (II) As(x) / 0,AG(x) / 
0, where in both cases Au(x) = Ad(x) = As(x) is assumed. Here the 
spin-dependent Altarelli-Parisi equations are applied directly to the gauge- 
invariant parton spin distributions. 

While in the GRSV sets, polarized deep inelastic lepton nucleon scatter- 
ing was considered up to the next-to-leading order QCD within the frame- 
work of the radiative parton model. The structure functions were subject 
to two different sets of theoretical constraints related to two different views 
concerning the flavor SUj(3) symmetry properties of hyperon /3-decays [16]. 
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One set is called the standard (ST) scenario, while the other set is called 
the valence (VA) scenario. 

Figure 2 shows the unpolarized differential cross sections being plot- 
ted with qx and \J~S = 500GeF. Three different dilepton masses of Q = 
WGeV, Q = 50GeV and Q = M z are used. The MRS(A') and the GRV 
sets, which are defined up to the next-to-leading order in QCD, give essen- 
tially the same results. 

Figures 3a-3c and 4a-4c show the predictions for the single-spin asym- 
metries and double-spin asymmetries, respectively. From Fig. 3a-3c, we 
note that the two CLW sets generally give larger differences between the 
single-spin asymmetries \Al\ than the two GRSV sets. The differences be- 
tween the asymmetries from the two GRSV sets are essentially constant over 
the given qx range. The asymmetries from both the CLW and GRSV sets 
are large when the dilepton mass is equal to Mz, due to the fact that the 
effects of the Z interference become important when the dilepton mass is 
equal to Mz- From Q = WGeV to Q = M z , the asymmetries change by 
two orders of magnitude. 

From Fig. 4a-4c, the double-spin asymmetries \All\ are of the same 
order of magnitude in the given q^ range and for the three Q values. The 
CLW sets give comparatively larger differences between the asymmetries in 
the same qx range than the GRSV sets, which give essentially the same 
differences between the asymmetries. 
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Finally, in Fig. 5a-5d, we give the contributions of the subprocesses G + 
q — > l + l~ + q (Fig. la-lb) and q+q — > l + l~ + G (Fig. lc-ld) to the single-spin 
and double-spin asymmetries at Q = Mz and = 500Gey. Fig. 5a-5b are 
obtained from the CLW sets, while Fig. 5c-5d are obtained from the GRSV 
sets. We observe that quark-gluon and quark-quark contributions to the 
single-spin and double-spin asymmetries behave differently with increasing 
qr- The contributions of the subprocesses to the asymmetries obtained from 
the CLW sets are more dispersed and are larger than those from the GRSV 
sets, while the single-spin or double-spin asymmetries of the subprocesses 
obtained from the two GRSV sets have similar behavior. In general, the 
contributions from quark-gluon (Fig. la-lb) vary more significantly than 
that from quarks (Fig lc-ld) alone. 

IV. Conclusion 

In this paper, we have made a numerical study of the single-spin and 
double-spin asymmetries for the PP — > l + l~ + jet process. In this study, 
we have used, for convenience, the helicity amplitude method to calculate 
the scattering amplitudes and the differential cross sections for the asym- 
metries. Both the asymmetries of the angular distribution of lepton pair in 
the large transverse momentum Drell-Yan process at RHIC energies were 
studied. Polarized parton distributions from CLW [12] and from GRSV 
[13] were used. Both the CLW and the GRSV sets were obtained to the 
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next-to-leading order. We have used the MRS (A') [14] and the GRV [15] 
sets of unpolarized distribution functions, which are also given to the next- 
to-leading order. The unpolarized differential cross sections obtained from 
these sets are essentially the same. 

We observe that the two CLW sets generally give larger differences be- 
tween the single-spin asymmetries \Al\ than the two GRSV sets, and the 
asymmetries generally increase as the dilepton mass approaches Mz, due to 
the large effects of the Z interference from the parity-violating Z-fermion 
coupling. The differences between the asymmetries from the two GRSV 
sets are essentially constant and small. From Q = WGeV to Q = Mz, the 
asymmetries change by two orders of magnitude. 

However, the double-spin asymmetries \All\ are of the same order of 
magnitude in the same qx range and for the three Q values. The CLW sets 
give comparatively larger differences between the asymmetries for the same 
qx range than the GRSV sets, which give essentially constant differences 
between the asymmetries. We note that the asymmetries are nearly constant 
for different Q values. 

We have also calculated the contributions of the subprocesses G + q — ► 
l + l~ + q (Fig. la-lb) and q + q -> l + l~ + G (Fig. lc-ld) to the single- 
spin and double-spin asymmetries at Q = Mz and y/S = 500GeV . The 
single-spin and double-spin asymmetries behave differently with increasing 
qx- The asymmetries of the subprocesses obtained from the CLW sets are 



15 



more dispersed than from the GRSV sets, while those obtained from the 
two GRSV sets have similar behavior. 

We observe that, at high energies, Z boson contributes significantly to 
the polarized process and induces large parity-violating single-spin asym- 
metries. The single-spin asymmetries become large as the dilepton mass 
approaches the Z peak for both the CLW and GRSV sets. The double-spin 
asymmetries are generally larger than the single-spin asymmetries for any 
dilepton mass. Therefore, measurements of the single-spin asymmetries may 
be used to pin down the spin structure of the proton for distinguishing the 
two CLW sets. However, the two GRSV sets give about the same results 
over large values of qr for different dilepton masses. 

Finally, we should mention that Leader and Srihdar [17] have also made 
a detailed study on the polarized PP — > l + l~ + jet and obtained some 
similar results. Their work and ours were based on the 0(a s ) tree Feynman 
diagrams on the parton level. One-loop radiative QCD corrections to the 
polarized parton subprocesses remain to be calculated and are difficult. We 
note that, for the unpolarized Drell-Yan process, inclusion of the 0{a 2 s ) 
corrections changes the cross section at fixed target energies up to 10% [18]. 
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Appendix A 



Here we list the nonvanishing helicity amplitudes for the polarized hadron 
subprocesses. With the contribution of a Z boson, the amplitudes read 
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The parameters L 9 = 2T^—2Q q sin 2 -R9 = — 2Q g sin 2 6*^/ are, respectively, 
the left-handed and right-handed coupling constants, of quark coupled with 
the Z boson, where T3 is an SU{2) isospin quantum number, Q q is the charge 



18 



of the quark, 9w is the Weinberg angle, and T a is the SU (3) generator with 
a being the color index. 

With the contribution of a virtual photon, we have 
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Appendix B 

In this appendix, we list the relations between the invariant variables Sij 
and the Mandelstam invariants. 
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S14 = 2pi ■ p A 



2 (Q " t) + 2(J^) COSa + T^ SmaCOS/? 



S15 = 2pi -p 5 = -U 
«23 = 2p 2 • P3 



= o <3 -« -TTF — 7^ cosa + ^ — -^5- sin a cos ^ 
2 2(s-Q 2 ) s-Q z 

s 2 4 = 2p 2 • P4 



l._ 2 QH-sii, ^QHtu . 

-(Q + — — -^cosa--^ — — sin a cos /? 

2 2{s-Q z ) s-Q z 



S25 = 2p 2 • P5 = -i 



S34 = 2p 3 • p 4 = Q 2 

u + t 

■S35 = 2p 3 -p 5 = — (1-cosa) 

u + i,^ 

s 45 = 2p 4 • p 5 = — (1 + cos a) . 

References 

[1] J. Ashman, et al. (EMC), Phys. Lett. B206 (1988) 364; Nucl. Phys. 
B328 (1989) 1. 

[2] D. Adams, et al. SMC, Phys. Lett. B302 (1993) 533; B329 (1994) 
399; B339(E) (1994) 332; Phys. Lett. B357 (1995) 248. 

[3] D.L. Anthony, et al. E142, Phys. Rev. Lett. 71 (1993) 959. 

[4] K. Abe, et al. E143, Phys. Rev. Lett. 74 (1995) 346; 75 (1995) 25. 



21 



[5] R.L. Jaffe and A. Manohar, Nucl. Phys. B337 (1990) 509. 

[6] G. Altarelli and G.G. Ross, Phys. Lett. B212 (1988) 883; R. Carlitz, 
J. Collins and A. Mueller, Phys. Lett. B214 (1988) 229; A. Efremov 
and O. Teryaev, Phys. Lett. B240 (1990) 373. 

[7] H.-Y. Cheng and C.-F. Wai, Phys. Rev. D46 (1992) 125. 

[8] H.-Y. Cheng, M. Huang and C.-F. Wai, Phys. Rev. D49 (1994) 1272. 

[9] P. De Causmaecker, R. Gastmans, W. Troost and T.T. Wu, Phys. Lett. 
B105 (1981) 215; D. Danckaert, P. De Causmaecker, R. Gastmans, W. 
Troost and T.T. Wu, Phys. Lett. B114 (1982) 203; P. De Causmaecker, 
R. Gastmans, W. Troost and T.T. Wu, Nucl. Phys. B206 (1982) 53; 
F.A. Berends, R. Kleiss, P. De Causmaecker, R. Gastmans, W. Troost 
and T.T. Wu, Nucl. Phys. B206 (1982) 61; B239 (1984) 382, 395. 

[10] R.D. Carlitz and R.S. Willey, Phys. Rev. D45 (1992) 2323. 

[11] G. Altarelli, R.K. Ellis and G. Martinelli, Nucl. Phys. B157 (1979) 
461. 

[12] H.-Y. Cheng, H. H. Liu and C.-Y. Wu, Phys. Rev. D53 (1996) 2380. 

[13] M. Gliick, E. Reya, M. Stratmann and W. Vogelsang, Phys. Lett. B359 
(1995) 201; Phys. Rev. D53 (1996) 4775. 



22 



[14] A.D. Martin, R.G. Roberts and W.J. Stirling, Phys. Rev. D50 (1994) 
6734; Phys. Lett. B354 (1995) 155. 

[15] M. Gliick, E. Reya and A. Vogt, Z. Phys. C67 (1995) 433. 

[16] H.J. Lipkin, Phys. Lett. B256 (1991) 284; B337 (1994) 157; 
J. Lichtenstadt and H.J. Lipkin, Phys. Lett. B353 (1995) 119. 

[17] E. Leader and K. Sridhar, Phys. Lett. B311 (1993) 324; Nucl. Phys. 
B419 (1994) 3. 

[18] P.J. Rijken and W.L. van Neerven, Phys. Rev. D51 (1995) 44, and 
references cited therein. 

FIGURE CAPTIONS 

Fig. 1 Feynman diagrams for the PP — > / + /~ + jet. Curly lines denote 
gluons, wavy lines denote 7* or Z bosons and solid lines denote quarks and 
lepton pairs. 

Fig. 2 Unpolarized differential cross sections for the PP — > l + l~ + jet 
process at as a function of qx, with r is fixed at Q 2 /S and 

y = 0, at various dilepton masses: Q = WGeV, 50GeV, M z . The MRS(A') 
set and the GRV set give essentially the same results. 
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Fig. 3 Single-spin asymmetries for the PP — > / + Z + jet process at 
as a function of qr- The asymmetries shown in these figures 
are constructed out of the cross-sections differential in qr, rapidity y and r, 
with r being fixed at Q 2 /S and y = 0. The full and dotted lines correspond 
respectively to the CLW Set I and Set II, which are used together with the 
MRS(A') unpolarized distribution functions. The circle-full and dotted-full 
lines correspond respectively to the GRSV standard (ST) and valence (VA) 
scenario sets, which are used together with the GRV sets. Three figures 
corresponding to (a) Q = lOGeV, (b) Q = 50GeV and (c) Q = Mz are 
shown. 

Fig. 4 Double-spin asymmetries for the PP — > l + l~ + jet process at 
= 500Gey, as a function of qr- Three figures corresponding to (a) 
Q = lOGeF, (b) Q = 50GeV and (c) Q = M z are shown. 

Fig. 5 Contributions of the subprocesses G + q — > l + l~ + q (Fig. la- 
lb) and q + q — > / + /~ + G (Fig. lc-ld) to the single-spin and double-spin 
asymmetries at Q = Mz and = 500Gey. (a) and (b) are contributions 
obtained from the CLW sets, while (c) and (d) are those obtained from the 
GRSV sets. 
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